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Abstract 



Qh. We provide the first explanation of the counter-intuitive scalar glueball cou- 

plings to pseudoscalar mesons found in lattice QCD and predict hitherto uncal- 
culated decay modes. Significant a\i: and (7T7r),s(7r7r)s couplings are found. We 
demonstrate the equivalence of linear and quadratic mass matrices for glueball- 
quarkonium mixing. The equivalence of formalisms which deal with a glueball- 
quarkonium basis and only a quarkonium basis is demonstrated. We show that 
the /o(1500) is not the heaviest state arising from glueball-quarkonium mixing 
for a glueball mass consistent with lattice QCD. The masses and couplings of 
scalar mesons, as well as their valence content, are calculated. 
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1 Introduction 



The existence of a gluon self-coupling in QCD suggests that, in addition to the 
conventional qq states, there may be non-gg mesons: bound states built from gluons, 
called glueballs. The abundance of qq meson states in the 1-2 GeV region and the 
possibility of glueball-quarkonium mixing makes the identification of the would-be 
lightest non-gg mesons extremely difficult. To date, no glueball state has been firmly 
established, although the existence of glueballs has been established in lattice QCD. 

Although the current situation with the identification of glueball states is rather 
complicated, some progress has been made recently in the scalar glueball sector, where 
both experimental and lattice QCD results seem to converge. Various lattice QCD 
glueball mass estimates have been made in the literature, and one of the differences 
stems from the way the physical results are obtained from the raw lattice data: either 
by calculating the sting tension, or the p mass. UKQCD estimates 1.55 ± 0.05 GeV 
by fixing to the string tension [TJ. GF11 originally estimated 1.74 ± 0.07 GeV M by 
fixing to the p mass. Later estimates include 1.71 ±0.06 GeV || and 1.65 ±0.06 GeV 
g. Attempts at reconciling UKQCD and GF11 results yielded 1.57 ± 0.09 @ and 
1.63 ± 0.09 ||. In what follows, we shall take the glueball mass tuq to be 1.6 GeV. 
Accordingly, there are two experimental candidates |§, / (1500) and /j(1710), in the 
right mass range. 

Recently, ref. JFj showed that the hypothesis where isoscalar meson mixing pro- 
ceeds through an intermediate glueball, called "glueball dominance" , can consistently 
account for isoscalar meson masses in various J PC sectors by employing glueball 
masses predicted by lattice QCD. Here we explore in detail the consequences of glue- 
ball dominance in the scalar sector, which differs from any other J PC sector due to 
the relative proximity of the glueball and quarkonia masses. Particularly, we demon- 
strate that the formulation of glueball dominance in refs. |7|, || is consistent with 
glueball-quarkonium mixing formulated in refs. H f|, |i~0| . Moreover, the recent 



lattice QCD calculation of ref. [[Ll|] found an unusual decay pattern for the scalar 
glueball, which is not consistent with any published model. We demonstrate that 
glueball dominance can explain this decay pattern. 

In Section 2 the canonical formulation of glueball dominance is introduced. Section 



2 



3 merges glueball dominance and the 3 Po/ flux-tube model to describe glueball decays. 
Section 4 discusses glueball-quarkonium mixing. Phenomenological implications are 
indicated in Section 5, and a summary given in Section 6. 



2 Glueball dominance 

We assume the glueball dominance of quarkonium mixing, viz., there is no direct 
quarkonium-quarkonium mixing, and the qq <-> qq transition is dominated by the 
glueball with the corresponding quantum numbers in the intermediate state [0 . 



Although the validity of glueball dominance has not been shown in QCD, the 
conclusions drawn from glueball dominance often agree with conclusions derived from 
the large number of colors N c limit of QCD: 

• Consider a Feynman graph where a qq pair annihilates into an arbitrary num- 
ber of gluons which then create a qq pair. In the large N c limit the graph is 
O(jj-) independent of the number of intermediate gluons. This corresponds 
to the finding (in glueball dominance) that the coupling between mesons via 
an intermediate glueball is largely independent of the C-parity of the glueball 
0, i.e. the number of gluons that a glueball can be built from in perturbative 
QCD. 

• Consider a Feynman graph where (a) two gluons each create a qq pair (i.e. two 
quark loops) which again combines into two gluons; (b) two gluons combine to 
a single intermediate quark loop and then combine to two gluons. (a) is 0{^) 
and (b) is O(jj-). Glueball dominance postulates that glueballs mix via a single 
intermediate meson, corresponding to (b), in agreement with the large N c limit. 

• The glueball dominance description of glueball decay postulates that the glue- 
ball couples to a meson, which then subsequently decays. In order for this 
process to happen, the time r for the qq pair created in the glueball to form a 
meson should be significantly less than the time 1/T required for yet another 
qq pair to form so that the meson decays. In the large N c limit, 1/T is 0(N C ). 
The time taken for the created qq pair in the glueball to form a meson should 
be inversely proportional to the mass of the state, so that r is 0(1). Hence the 
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large N c limit gives the condition r < 1/r required by the glueball dominance. 
Another way to see the large N c result is to note that quark pair creation is sup- 
pressed (O(jj-)) but quark rearrangement is not ((9(1)), so that created quarks 
in a glueball would rearrange to form a meson, rather than create another qq 
pair. 

Here we review the glueball dominance picture discussed in more detail in ref. 
0. The possibility of the transition qq — > gg . . . g — > qq is accounted for by the 
quark mixing amplitudes, A qq i, which are included in the meson mass square matrix 
(written down here in the ss, (uu + dd)/y/2 basis): 



m 2 sS + A ss V2A sn 
V2A ns + 2A 



(2.1) 



where m sS and m nn , n = u,d are the primitive (bare) quarkonia masses. 
The quark mixing amplitudes can be represented in the standard form, 

A ~ (Vertexl) ■ (Propagator) • (Vertex2), 

which reduces to 

(qq\Hll\i)(i\H«i]q'q') 
Y M 2 (q'q) -M 2 (z) ' 1 ' 

where H^ c is the quark pair creation operator for the flavor q, and \i) is a complete 

set of the (gluon) intermediate states. Because of the assumed glueball dominance, 

the sum (2.2) is saturated by the low-lying glueball: 

A qq ~ a , f gqG = (qq\H^\G) 2 , (2.3) 

for q = n(= u,d),s, and, in view of the factorization hypothesis discussed in more 
detail in ref. 0, 

A sn = A ns = \j A nn * A ss ~ , (2.4) 

\JVKn- m G ){ m ss- m G ) 

where f qq -G defined in (2.3) is the coupling of the intermediate glueball to qq, and m G 
is the corresponding (physical) glueball mass. 
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Thus, the mass matrix (2.1) reduces to 



f 2 

2 _|_ JssG 



VZfnnGfssG 



,2 



m nn + mi. -mi 



(2.5) 



/ 



The masses of the physical isoscalar states fo and / are obtained by diagonalizing 
this mass matrix: 



( m 2 , 

Jo 



V 



o 



(2.6) 



m 



fo 



3 Glueball decay 



3.1 Scalar glueball decay to two pseudoscalar mesons 

Glueball decay via glueball dominance was introduced in ref. ||13| . We follow their 
approach, except for one improvement. Ref. |]TB| coupled the decaying glueball to 
an off-shell meson, which then subsequently decays to the outgoing mesons. The 
coupling used for this latter process is extracted from experiment where the off-shell 
meson is on-shell. This should not introduce undue errors unless the coupling is 
strongly dependent on energy of the off-shell meson. 

In this section we shall deal with a case where the coupling is strongly dependent 
on the energy of the off-shell meson. Our proposed solution is to calculate the decay 
for the correct energy of the off-shell meson by employing the phenomenologically 
successful 3 Pq model [14, lq| . Since the nonrelativistic 3 Po and Isgur-Paton flux- 
tube decay model give identical predictions for simple harmonic oscillator meson wave 
functions |16|], which we employ, our predictions can also be viewed as predictions of 
the flux-tube model. 

In fact, the flux-tube model, motivated from the strong coupling limit of the 
Hamiltonian formulation of the lattice gauge theory (HLGT), added to glueball dom- 
inance affords an intuitive picture of the decay process. In the lowest order in pertur- 
bation theory glueballs can be viewed as rings of flux in HLGT. Pair creation occurs 
in the first order of perturbation theory and breaks the flux-ring up into a meson. 
The flux-tube connecting the two quarks in the meson then breaks via the creation 
of a quark-antiquark pair with vacuum ( 3 Pq) quantum numbers to form two outgoing 
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mesons |Tj 

The amplitude (in GeV) for the decay of the scalar glueball to two outgoing 
mesons is then given by 

■ j V^fnnG a 1 -\ , fssG a f -\ (n 1 N 

M = T A{nn) H T A{ss), (3.1) 

m nn ~ m G m ss ~ ™ G 

where f nn o and f sS G are the couplings of the scalar glueball to the intermediate nn 
and ss scalar mesons, respectively, introduced in Section 2. We have also taken 
care to insert the scalar meson propagator with the masses ordered according to the 
prescription of glueball dominance (2.5). A(nn) and A(ss) are the intermediate scalar 
meson 3 P model decay amplitudes to two outgoing mesons. We have assumed that 
only the ground state scalar mesons saturate the decay. 

The full 3 Po model amplitude is given in Appendix B of ref. |14]]; here we just write 
down the case of identical inverse radii Pa = Pb — Pc = P f° r identical quark masses 1 , 
for simplicity (A denotes the scalar meson, and B and C the outgoing mesons): 

| 8m 2 G M B M c lforl 2p 2 p 2 

A = i m a 70 (1 - w ] exp{ ~i2^ } ' (3 - 2) 

where we neglected the factor arising from the flavors of the mesons. The mock 



meson phase space convention is specified by Mj^b,c 01) an d p is the momentum 
of the outgoing meson B in the glueball rest frame. The width is computed from 
the amplitude in Eq. ( |3.2j ) by using the standard formula given by the Particle Data 
Group [0. The pair creation constant 70 is usually taken to be the same for nn and 
ss pair creation |L4|, [TEfl , as we shall do here. 
The composition of n and rj : is 

rj = sin 6 \nn) + cos 6 \ss), n = cos 8 \nn) — sin 8 \ss). (3.3) 

We obtain the following simple relationships between the 3 Po model amplitudes 

when Pa, Pb, Pc,P and the mock meson phase space are taken to be constant for all 

meson decay processes: 

,,uu + dd , t/ uu + dd T ^ T ->, ..uu + dd . 

► tttt) : — > KK) : A{— -y= > nn) : 



-'^The 3 Po model amplitudes depend explicitly on the light and strange quark masses for decays 
where the initial quarks are different from the quarks in the created pair. In this work we take the 
light and strange quark masses to be identical for decays to KK. 
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A( UU "t-^ > r]r] ) : A(ss — ► fTTT) : ^4(ss — ► 7777) : ^4(ss — ► 7777 ) 

V2 

= 2 : 1 : 2 sin # : 2 sin cos 6 : V2 : 2 v / 2cos 2 # : -2^2 sin 6* cos 0. (3.4) 
When we take m /""^ = m /^y (explained below), we obtain 



.M(G^7r7r) : TW(G^KiT) : M{G ^ r)rf) : M{G ^ r]r}') = 

1:1:1:0. (3.5) 

This is the result one obtains when naively coupling the quarks in the outgoing 



mesons to the vacuum [TjJ], often referred to as "flavour democratic coupling". It 
would also yield a horizontal line for our predicted amplitude in Fig. 1. 

The lattice results were obtained in the SU(3) limit. To compare we shall also 
adopt the SU(3) limit in the remainder of this section. Hence we take the couplings 
and quark masses to be identical, i.e. f nfl G = fssG = fsu(3)G and m nfl = m sS = m S u(3)- 

Fig. 1 shows our results. 

The solid line represents our basic prediction. We use (3 = 0.4 GeV found to 
enable a fit of a large range of meson decays |nj [T8|| . Mock meson phase space is 
employed since this enables a prediction of M. for all p, as can be done in lattice 
QCD. Since we work in the limit of 577(3) symmetry, we take the mock meson phase 
space parameters to be those of say ss mesons, i.e. Ma = 1.49 GeV, Mb, Mq = 0.85 
GeV (see Table 5 of ref. fl4]l ). We take the pair creation constant to be 70 = 0.39 

0- 

Since all parameters are constrained, except for fsu{3)G/( m su(3) ~ m <§)' we re g ar d 
our fit in Fig. 1 as a one-parameter fit. This parameter sets the overall decay 
strength of the glueball, and does not influence the relative strengths of the various 
decay modes. 

The lattice results^ for the 7777, KK and nir decay modes are plotted from left to 
right as the data points. It is non-trivial that our results are consistent. The fact that 
the prediction in Fig. 1 is not a horizontal line, as one naively expects, indicates that 



2 GF11 predicts ^ = 0.834t^?|, 2.654±g;^, 3.099±g;^g for scalar glueball decay to tttt, KK 
and 7777, respectively 11, ^o), where m p is the mass of the p meson. The predictions are given as a 



function of pseudoscalar mass, which we translate to the momentum of the outgoing mesons p using 
conservation of energy: p 2 = wlq/A — m 2 PS , where mps is the relevant pseudoscalar meson mass. 
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the detailed dynamics of the flux-tube and 3 Po models combined with the hypothesis 
of glueball dominance captures the correct strong interaction dynamics. This success 
is not shared by other models of glueball decay based on perturbative QCD decay 



dynamics, where the naive pattern of Eq. ( |3.5|) arises |21| . The other points and lines 
in Fig. 1 indicate parameter variations and are discussed in the caption of the figure. 
We fit 



fsU(3)G 



2 2 
m G ~ m SU{3) 



0.34 ±0.04. (3.6) 



3.2 Scalar glueball decay to two mesons 

Having predicted the scalar glueball decay to pseudoscalar mesons, we are now in 
a position to make the first predictions in the literature of the decay of the scalar 
glueball to non-pseudoscalar mesons. 

Since we have fitted Eq. (3.6) using mock meson phase space, we again use this 
convention and hence 70 = 0.39 JT4|. We again use (3 = 0.4 GeV and do the calculation 
in the SU(3) limit with the ss mock meson masses^ [|14]]. The results are indicated in 
Table 1. The primitive glueball amplitudes should be understood to be correct up to 
a sign. The analytical expressions used for the amplitudes can be found in Appendix 
A. 

We see from Table 1 that the total width of the 1.6 GeV scalar glueball is 250 — 390 
MeV excluding (im) s (mi) s decays. There is also substantial phase space dependence 
for the glueball decay amplitudes to tttt, KK, rjr] and 7r(1300)7r. 



4 Glueball mixing 

4.1 Glueball— quarkonium basis 

3 Taking the parameters for nn mesons, i.e. Ma = 1-25 GeV, Mb : Mc = 0.85 GeV would 
give scalar meson widths to irn, KK and 7777 ~ 20% larger. A fit to the lattice data then yields a 
different value for /sc/(3)G/( m |(7(3) _ m ^)- The predictions for the glueball widths to irir,KK and 
7777 are identical, and the dominant width in Table 1, to ai7r, is 1% different from the value quoted. 
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4.1.1 Glueball— quarkonium linear mass matrix 

In ref . [[KJ , Weingarten suggested the following 3x3 linear mass matrix, which stems 
from the Hamiltonian formulation of QCD, to describe the mixing of a glueball and 
quarkonia: 

V2z \ 

z m sS , (4.1) 
V2z 

Winn J 

where z stands for the annihilation amplitude of quarkonium into a glueball which 
has dimensionality (mass) and represents a counterpart of our /'s which have dimen- 
sionality (mass) 2 . In order to test our results by comparing with available lattice QCD 
data, we should establish a relation between this linear mass matrix and our mass 
squared one. 



V 



4.1.2 Glueball— quarkonium quadratic mass matrix 

We first rewrite Weingarten's matrix for the squares of the glueball and quarkonia 
masses and show its equivalence to glueball dominance in the 2x2 subspace spanned 
by quarkonia. We then establish a relation between Weingarten's linear and our 
quadratic mass matrices. 



So, consider 



rn 



a 



f y/2f 



f 
V2f 



mz 







mi 



\ 



J 



(4.2) 



where the vanishing off-diagonal elements indicate that there is no direct quarkonium- 
quarkonium mixing, i.e., glueball dominance. 

Proposition 1. The mass matrix (4.2) is equivalent (gives the same physical 
quarkonia masses) to glueball dominance in the 2x2 subspace spanned by quarkonia 
(with f nnG = f sSG ). 

Proof. First, we rewrite the mass matrix (2.5) (with f nn G — fssc) i n the following 
form: 

/ 



m ss ~ 7,2 A 



\ 




A = 



f 1 
j 1 



nnG 



m G - m ln 



r = 



\ 



m G ~ m ln 



m 



G 



mi 



(4.3) 
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where r is a complex number. The masses of the two physical states are now deter- 
mined from the equations 

m%+m 2 f , = m 2 nn + m 2 sS -A(2 + r 2 ), (4.4) 

m \ m \ = m nn m L- - A ( 2m L- + r2m ln) • ( 4 - 5 ) 

We take the equivalence of the matrices (2.5) and (4.2) to mean the equality of 
the corresponding eigenvalues m 2 f , m 2 ,. 

J0 Jo 

The eigenvalues of (4.2) are determined from the following three equations: 

m % + m /o + m \ = m G + m ln + mis, ( 4 -6) 

22,22,22 22,22,2 2 of2 /< m 

m G m fo + m G m f + m fo m f' = m G m nn + m G m ss + m nn™ss ~ 3/ , (4.7) 

m2 G m )o m \ = m G m ln m2 ss ~ f ( 2m L + m «n) , ( 4 - 8 ) 

It then follows from Eqs. (4.4), (4.6) and Eqs. (4.5), (4.8) that: 

_ 2 (2m 2 sS + r 2 m 2 n 2 (2 + r 2 )m^m^\ , , A n . 

f = O 2 7 2 K ~ o 2 T 2 A > (4- 9 ) 



and from Eqs. (4.4)-(4.7): 

2m^ + r 2 m 2 
3 " u 3 



/2 = | ^_^ m | _ ^nnT> ^ j A (4 1Q) 



Using the definition of r in Eq. (4.3), the equivalence of Eqs. (4.9) and (4.10) follows 
by simple algebra. Also, using the definition of A in Eq. (4.3), and inserting it in 
either Eqs. (4.9) or (4.10), we obtain that ±/ = f n nG — fssG- ^ 

This important result means that glueball dominance is nothing else but an effec- 
tive representation of the glueball-quarkonia mixing in the 2x2 subspace spanned 
by quarkonia. The relation is only possible because both formulations descibe the 
physics in terms of the same basis states. 

It is natural to define / = f nfl G = fssc , noting that Eq. (4.2) assumes 577(3) 
symmetric couplings. This definition is consistent with what is obtained when the 
equivalence of the 2x2 and 3x3 formalisms is demanded in Proposition 1. 

The mass matrix (4.2) possesses, however, more generality than the naive glueball 
dominance picture in the 2x2 quarkonia subspace. This is because the former, in 
contrast to the latter, allows one to obtain the valence glue content of the physical 
quarkonia, and the valence content of the physical glueball. 
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Proposition 2. The linear, (4.1), and quadratic, (4.2), formulations for the scalar 
mesons are equivalent provided that (i) z 2 <C m G , m 2 ^, m 2 s , (ii) m S s—m nn <C m nfi , m s5 
or m nn , m sS < m G . 

Proof. We take the equivalence of the linear, (4.1), and quadratic, (4.2), formula- 
tions to mean (i) the equality of the eigenvalues of the matrix (4.2) to the eigenvalues 
squared of the matrix (4.1), (ii) the equality of the eigenvectors of both matrices, for 
the same values of the input parameters mc, fn nn , m sS . Denote the matrices (4.1), (4. 2) 
by Mii n and M qaud , respectively, the corresponding diagonalized matrices by and 
A-qaud, and the matrix that diagonalizes the M's by S (it is the same for Mi in and 
M qau d because both have by construction the same eigenvectors). 

Since A Un = diag (m G ,m f >,m fo ), 

J 

A L = dia g ( m G , m2 f ', m %) = \aud- 

J 

It follows from this relation and 

^■un = SMi in S 1 ■ SM Hn S 1 = SMf in S 1 , A qaud = SM qaud S 1 

that 

M qaud = M? m . (4.11) 

Thus, the linear and quadratic formulations are equivalent provided that Eq. (4.11) 
is valid. 

Since the square of the mass matrix (4.1) is 



1 m 2 G + 3z 2 z(m G + m sS ) ^/2z(m G + m nn ) X 



z(m G + m sS ) m 2 sS + z 2 \f2z 2 
V2z(m G + m nn ) V2z 2 m 2 nn + 2z 2 y 



(4.12) 



it is clear that Eq. (4.11) follows if both conditions (i) and (ii) of Proposition 2 are 
satisfied. □ 

The equivalence of (4.2) and (4.12) also implies the following relation between z 
and f:f — z{m G + m qq -). In the scalar sector where the glueball and quarkonia have 
comparable masses, it reduces to 

/ ~ 2zm G . (4.13) 
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For any other J PC multiplet, mc 3> m q q, and, respectively, / ~ zuiq- 

An explicit numerical example on the equivalence of linear and quadratic mass 
matrix formulations is given in Appendix B. 

4.2 Scalar meson spectroscopy 

As the relation between Weingarten's linear mass matrix and our quadratic mass 
matrix is established in the previous subsection, we are ready to consider scalar 
meson spectroscopy implied by glueball dominance, and compare our results with the 
lattice QCD simulations of refs. [^, [Tt| . 



We shall first show that within the glueball dominance hypothesis, the /o(1500) 
cannot be the heaviest isoscalar scalar meson arising from ground state nn, ss and 
glueball mixing if mg > 1.5 GeV.f] 

The argument is as follows. There are three possibilities: (i) m sS < ma, (ii) 
m sS > mo, (iii) m sS = ttiq. The main property of the 3x3 mass matrices (4.1), (4. 2) 
(we do not prove it here) is that upon mixing the higher mass primitive state becomes 
more massive, while the lower mass primitive state becomes less massive (i.e., the 
mass splitting between the higher and lower mass primitive states increases as a 
result of the mixing). Therefore, in the case (i) uiq > nic > 1.5 GeV; in the case 
(ii) m f ' > m sS > mc > 1.5 GeV. Finally, in the case (iii) it can be shown that the 



physical ss and glueball states have masses ~ \j m G ^ /> anc ^ therefore, one of them 
is always higher than mc > 1.5 GeV. 

Hence, the /o(1500) is not the heaviest isoscalar scalar meson arising from ground 
state nn, ss and glueball mixing. If both the existence of / (980) and / (1370)// (400— 
1200) are confirmed by experiment, there has to be an extra degree of freedom to 
account for the existence of these states. 



4 This conclusion may be modified by the inclusion of effects beyond glueball dominance, e.g., 
coupled channels |2!q] , which may be especially relevant for scalar states [^3|. Ref. |2^] finds that 
the masses of the states are always lower than those of the primitive states. We believe this to be 
an artifact of the inclusion of only low-lying channels, which could lead to misleading results p3] . 
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4.2.1 The glueball— quarkonium coupling / 

In ref. ||, Lee and Weingarten estimate the mixing parameter z to beQ 

\z\ = 56±37MeV. (4.14) 

Eq. (4.14) implies, via (4.13) with m G = 1.65 ± 0.06 GeV §, 

\f\ ~ 0.19 ±0.13 GeV 2 . (4.15) 

This implies that in the flavor SU(3) limit where we take msu(3) — m nn — m ss — 
t^sn = m K* = 1-43 GeV || (and the same m G = 1.65 ± 0.06 GeV), 

/ 



2 2 
m G — m SU(3) 



(4.16) 



consistent with the value needed from glueball decay (Eq. (3.6)). 
We take from Eq. (3.6), defining m q q = {m nn + m sS )/2, 

/ 



0.34 ±0.04. (4.17) 



4.2.2 Two simulations for glueball— quarkonium mixing 

We now wish to consider two simulations for the glueball-quarkonia mixing based on 
the quadratic mass matrix (4.2), and extract the masses of the primitive quarkonia 
and the physical states, with the help of Eqs. (4.15) and (4.17). 

We use Eqs. (4.6)-(4.8) for the masses of the physical states m G , m/ , rrij>. 

For the first simulation, we also employ Eq. (4.17) and 

m G = 1.6GeV, m sS - m nn = 250 ± 50 MeV, (4.18) 

the latter being a typical mass splitting^] between ss and nn states for different meson 
5 In a most recent paper [Q, Lee and Weingarten introduce SU{2>) breaking effects, in terms of 
different values of z for G-nfi, G-ss mixing, ZGnh/zGss = 1.198±0.072, and estimate \za s §\ = 43±31 
MeV, which implies, with the above ratio, |zGnn| = 54 ± 40 MeV, similar to (4.14). 

6 Note that we could use the mass squared splitting for quarkonia, m 2 sS — to^ s , in place of the 

linear one in Eq. (4.18). With, e.g., the value for this mass squared splitting 0.65±0.01 GeV 2 , since 
for the remaining three P-wave nonets (in GeV 2 ) 2(m 2 K , — m 2 2 ) ~ 0.64, 2(m 2 KiA ~ m 2 ^) ~ 0.66, 



m r , 



2(m 2 KiB - m 2 bi ) ~ 0.65 gj, the solution is m d = 1694 ± 27 MeV, m fo = 1242 ± 21 MeV, 
1314 ± 12 MeV, m sS = 1542 ± 8 MeV, |/| = 0.285 ± 0.06 GeV 2 , and ^{m 2 nfl + mj g )/2 = 1432 ± 10 
MeV, consistent with the case of the linear mass splitting (middle column, Table 2). 
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multiplets || p^j , and solve the system of 5 equations (4.6)-(4.8),(4.17),(4.18) for m^, 
m/ , m nfi , m sS and /, by fixing m f > = 1.5 GeV. 

The reason for the latter requirement is that /o(1500) is established in more decay 
channels that any other scalar meson, and we should therefore construct our simula- 
tion of scalar meson spectroscopy with the constraint that one of the masses of the 
physical states is 1503 ± 11 MeV f§. 

For the second simulation, we employ Eq. (4.15) in place of (4.17). The solution 
to Eqs. (4.6)-(4.8), (4.17), (4.18) for the first case, and Eqs. (4.6)-(4.8), (4.15), (4.18) 
for the second case are presented in Table 2.[] 





Eqs. (4.6)-(4.8), (4.17), (4.18) 


Eqs. (4.6)-(4.8), (4.15), (4.18) 


m d , MeV 


1703 ±40 


1649tg 


m /o , MeV 


1218 ± 70 


1248l|| 


m sS , MeV 


1546 ± 17 


1527 ±23 


m nn , MeV 


1296 ± 33 


1277 ±27 


|/|, GeV 2 


0.305 ±0.09 


0.19 ±0.13 



Table 2: Solution to Eqs. (4.6)-(4.8), (4.17), (4.18), and Eqs. (4.6)-(4.8), (4.15), (4.18). 



We note that the value of / obtained in the first case is consistent with the value 
extracted from lattice QCD (Eq. (4.15)). It is also in agreement with values extracted 
phenomenologically for different J PC meson nonets in ref. [7|] which all are in the 
interval 0.27 — 0.32 GeV 2 . We however disagree with lattice QCD that the primitive 
ss is at least 200 MeV below the primitive glueball @, g but only ~ 70 ± 30 MeV, 
as seen from our solutions for m sS in the two cases considered. We also note that 
\J ( m nn ± m ss)/2 which is 1426 ±24 and 1407 ±25 MeV in the two cases, respectively, 
is consistent with rriK* = 1429 ± 6 MeV || . 

7 Note that, although we do find a solution with m sS > tuq for both simulations, the values 
for m„n and m sS obtained are too high to be accommodated by any of the existing quark models; 
typically m n n ~ 1.6 GeV, m sS ~ 1.8 GeV. Also, /o(1500) is the lightest of the three scalars. Our 
results are therefore in agreement with the conclusion of ref. || that the situation where the primitive 
ss state has a higher mass than the primitive glueball is incompatible with lattice QCD. 
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With / > 0, the valence content of the three physical states obtained in the first 
simulation, is 

1 1703) = (0.821±0.02)|G) + (0.493±0.02)|gg) + (0.287±0.05)| MM ^ rf ), (4.19) 

v2 

|1500) = -(0.410±0.04)|G) + (0.860±0.02)|ss)-(0.305±0.08)|^^^), (4.20) 

V2 

|1218) = -(0.397±0.08)|G) + (0.133±0.05)|ss) + (0.908±0.05)|^^^), (4.21) 

V2 

and shows that the physical glueball contains ~ 70% glue and ~ 30% qq, while each 
of the physical quarkonia contains ~ 15% glue and ~ 85% qq. The overall signs for 
the states have no physical significance. 

Although the masses of the physical states do not depend on the sign of /, the 
valence content of the physical states does. Namely, we find that under the inversion 
of the sign of / both the quark content of the physical glueball and the glue content 
of the physical quarkonia change their sign. However, it is not difficult to see that the 
3 Pq model decay width of the physical states remains invariant under the inversion of 
the sign of /. Hence, for the study of masses and decays of scalar mesons, one need 
to consider the case / > only. 

Notwithstanding the similarity of the results obtained, there is a principal differ- 
ence between our approach and that of refs. || f|, [lOfl . Lee and Weingarten choose 
the input parameters uig, rn s g, m n n an d z to obtain the three physical masses. Al- 
though the input parameters vrtQ and z obtained are consistent with their lattice QCD 
calculations, the mass splitting m sS — m nii ~ 50 MeV, is atypical of that provided 
by constituent quark models in Eq. (4.18). Moreover, the Gell-Mann-Okubo type 
relation m^ + m^ = 2m 2 sn = 2m 2 K * [^4[ does not hold in their approach. In contrast, 
our approach does not assume m uU and m sS individually, but only their difference. To 
compensate, we need one more input parameter: we require that one of the physical 
masses be in agreement with the well established scalar state /o(1500). 

5 Phenomenology 

For some time there has been a controversy over the spin J assignment of /j(1710), 



and hence the existence of a scalar state at this mass |19]. This matter is not yet fully 
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resolved. Our results in Section 4.2 suggest that /o(1500) cannot be the heaviest state 
arising from glueball-quarkonium mixing. Hence we highlight recent evidence for a 
J = component in /j(1710). BES seperated both a J = and a J = 2 component, 
with the scalar state having mass 1780 MeV and a width of 85 ± 25 MeV ||25| . There 
are also claims of a J = signal at 1750 ± 15 MeV with width 160 ± 40 MeV PB 



and with mass 17041^ MeV with width 124±fl MeV f27j. m & obtained in Table 2 is 
consistent with these experimental masses. 

The ai7r decay of the primitive glueball is expected to be larger than any pseu- 
doscalar decay mode. There is some evidence for the production of / (1500) via af 
exchange in the reaction n~p — > 7T + 7T~n, i.e. an a\"K coupling of the /o(1500), in 
CERN-Cracow-Munich data with a polarized target |28|| . We urge experimenters 
to allow for the a^n decay mode in partial wave analyses. This applies to Crystal 
Barrel at CERN for analysis of 7r + 7r~7r°7r° [^J and 7r + 7r~7r + 7r~ data, and to Mark III 
7r + 7r~7r + 7r~ data p6| . 

The (7r7r)s(7r7r),s decay of the glueball can be very substantial, depending sensi- 
tively on the width and mass of the intermediate fo coupling to (t[tt)s (see Table 1). 
There are indeed indications of substantial (7T7r)s(7r7r)s decay modes in / (1500) and 
/o(1710) |§§. 

The two-photon widths of the various states provide stringent consistency checks 
for our results. In the flavour SU(3) limit the 77 width for a state ac\G) + + 



a,, 



uu +^ d ) i s proportional to 
x " 2 



v s5 + "7f an ™J ^ 5 '^ 

where we have taken the charges of the quarks into account, and normalized the 
expression to be unity when the state is built purely from ss. 

For the states in Eqs. (4.19)-(4.21) we obtain the two-photon width ratio for the 
states /o ( 1 370), /o ( 1 500), /o( 1710) to be 11 : 0.05 : 2. 

If one takes model-dependent estimates of two-photon widths of ss, i.e. 0.16 - 
0.20 keV fn| , one observes that the two-photon width of /o(1500) is consistent with 
the experimental bound of < 0.17 keV f32fl . Since /o(1370) is dominated by light 
quarks, our estimate for the two-photon width of /o(1370) is ^ 11 times the two- 
photon width of (3.25-6.46 keV [§TJ). This is consistent 5.4 ±2.3 keV § from 
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experiment. There are currently no experimental estimates for 77 width of the J = 
component of /j(1710). 

From Eqs. (4.19)-(4.21) it is clear that the glueball content of /j(1710) (J = 0) 
is significantly higher than that of /o(1500) and /o(1370), which are similar. This 
can be tested by evaluating the states' coupling to two gluons in the model of ref. 
19|: B R(f q(1710) -> gg) > 0.5 Q and Si2(/ O (1500) -> gg) = 0.3 - 0.5 or 



0.64 ±0.11 m. Since the expectation for a glueball is that BR(G -> ##) > 0.5 |4 



both / (1710) and /o(1500) are consistent with a sizable glueball component, and 
with / (1710) having a larger glueball component. Moreover, ref. |35| concludes that 
/o(1370) may have some glueball admixture, smaller than / (1500) and /o(1710), but 
is dominantly quarkonium, partially in agreement with our results. Predictions here 
are complicated by the large width of / (1370) fl3~4f . 



6 Summary 

In this paper we suggest a coherent view at the scalar glueball as having the following 
properties: 

(i) A (physical) intermediate state in scalar qq annihilation-creation transitions 
(called "glueball dominance"). 

(ii) A state decaying to two mesons via an intermediate scalar meson. 

(iii) A primitive state which mixes with the primitive nn (n = u, d) and ss quarko- 
nia to form three physical scalar mesons. 

Three main assumptions are employed in this work: glueball dominance, SU(3) 
symmetry and the assumption that only ground state quarkonia are relevant to scalar 
glueball mixing and decay. 

As can be seen in Figure 1, glueball dominance together with the calculation of 
energy dependent couplings in the 3 Po/flux-tube model can account for the counter- 
intuitive primitive glueball couplings to tctc, KK and rjrj found in lattice QCD. In 
Table 1 a total glueball width of greater than 250 — 390 MeV with a dominant a^Tc 
decay of 70 — 180 MeV is predicted. Decay to (7r7r),s(7r7r),s may also be significant as 
was observed experimentally for / (1500) and /j(1710). 

The quadratic mass matrices in the 2x2 quarkonium and 3x3 glueball-quarkonium 
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formulations are equivalent. The linear and quadratic 3x3 glueball-quarkonium mass 
matrices are equivalent under the requirements that (i) z 2 <C m^m^m^ and (ii) 
m sS — m n n <C m nf i,m S s or m n n,m sS <C ma- The conditions mentioned are always 
fulfilled in this work, and an illustrative example of the equivalence can be found in 
Appendix B. 

The /o(1500) is not the heaviest state arizing from glueball-quarkonium mixing, 
implying that if the existence of both / (980) and /„(1370)// (400 - 1200) is con- 
firmed, allowance should be made for an additional degree of freedom. The glueball- 
quarkonium coupling extracted from our glueball decay analysis is consistent with 
estimates from lattice QCD. 

The mass of the physical glueball is consistent with the experimental /j(1710). 
Experimental two-photon and J/ip radiative decay data are consistent with the va- 
lence content predicted for the physical states. 
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Appendix A: Glueball decay couplings 

The glueball decay amplitudes (evaluated in Table 1) to various outgoing states are 
(for Pa = Pb = Pc = P and identical quark masses) 



M = 




I \ M A V P 



8m 2 G M B M c [2 



7T4 7 exp{ 




} G7, 



m 



5(7(3) 



where 




for pp, ojijj, K*K* 



for 7T7T, KK, 7]7] 

for pp, uu, K*K* 



(D-wave) 



(S-wave) 
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zu = -^-7! for a\TT (P-wave) 

w = |- (l + ^) 2 - ^) 4 ) for ^(1300)^ (S-wave) 

*> = ™ (l + ^(|) 2 " ^(|) 4 ) for /o(600)/ (600) (S-wave) 



Appendix B: Numerical example of the equivalence 
of linear and quadratic mass matrix formulations 

An illustrative example of how Proposition 2 works is the case analyzed by Weingarten 



10 1 where for the input masses m nn = 1450 MeV, m s5 = 1516 MeV, mg = 1635 MeV 



and z = 77 MeV, the masses of the three physical states are tuq = 1710 MeV, 
77V — 1-5 GeV and mr n = 1390 MeV. To translate this linear mass case into our 

Jo J 

mass squared one we use the relation / = z(tog + m qq -) near Eq. (4.16) which, with 
m qq = (m n n + m S g)/2, gives / ~ 0.24 GeV 2 . When this / and the same input masses 
squared are used in (4.2), we obtain the masses of the physical states, rriQ = 1706 
MeV, to,' = 1.5 GeV, to/ = 1383 MeV, which are in excellent agreement with 
Weingarten's case, and the valence content of the physical states, 

, „. , ,uu + dd, 

|1706) = 0.87|G) + 0.34|ss) + 0.36| ^-), 

v2 

|1500) = -0.18|G) + 0.89lss) - 0.41| UU + 

v2 

, „. , , uu + dd, 

|1383) = -0.46|G) + 0.29|ss) + 0.84| ^-), 

v2 

which is in excellent agreement with the corresponding valence content of ref . jnj : 



|1710) = 0.871 G) + 0-34|gg) + 0.36| + J' d ), 



V2 

|1500) = -0.19|G) + 0.90|gg) -0.40| MM ^ ), 

v2 

, „. , .uu + dd, 

|1390) = -0.46|G) + 0.28|ss) + 0.84| -^-). 

v2 
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Fig. 1: The amplitude A4 (in GeV) plotted against the square of the momentum in the 
outgoing state p 2 (in GeV 2 ). The sold line is our basic prediction and the data points are 
the lattice predictions, both of which are discussed in the text. We allow the following 
parameter variations, (a) The dashed line differs from the basic prediction in that we take 
inverse radii motivated from studies of masses and wave functions of mesons and glueballs, 
not decays. The ir and K inverse radii were estimated as 0.54 GeV and 0.53 GeV respectively 
[j36"f , so we take (3b = (3c = 0.54 GeV. We also talk the glueball to have an r.m.s. radius 
\J (r 2 ) of ~ I fm [3?]] . Assuming that the scalar (P-wave) meson coupling to the glueball 
has the same size, one estimates for S.H.O. wave functions that (3 Ay/ ( r ' 2 ) = ^/§> yielding 
Pa = 0.6 GeV. (b) The solid black dots differ from the basic prediction, in that we adopt 
the relativistic phase space convention |l8|, Here we make the replacements Ma — > tug 
and Mb,Mq — > \Jvn PS + p 2 in Eq. (|3.2|) , where nips is the outgoing pseudoscalar mass. 
From left to right the points correspond to the nn, KK and tttt decay modes. The large 
points correspond to 70 = 0.4 |18| , and the small points to 70 = 0.53 jig] . 
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M 




Decay Mode 


Wave 


1.4 GeV 


1 r /~1 t r 

1.5 Gev 


l.o GeV 


1.7 GeV 


1 O ,,\ T 

l.o GeV 


Width (MeVJ 


7171 


c 



1.31 


0.96 


0.55 


0.093 


-0.41 


6 (6) 


KK 


S 


2.8 


2.7 


2.0 


1.6 


1.0 


81 (82) 


7777 


s 


3.2 


2.9 


2.4 


2.0 


1.4 


27 (27) 


a\7i 


T ) 

r 


10. 7p 


11. 2p 


11. 7p 


12. Op 


12. 3p 


177 (67) 


7T(1300)7T 


S 




1.1 


1.4 


1.8 


2.3 


23 (7) 


PP 


s 






3.3 


3.0 


2.6 


61 (46) 


ujijj 


s 






3.5 


3.1 


2.7 


16 (12) 


K*K* 


s 










4.1 




(7r7r) 5 (7T7r)s 


s 


< 9.6 


< 10.2 


< 10.9 


< 11.5 


<12.1 


< 490 (< 160) 



Table 1: Amplitudes for the decay of a scalar glueball to two mesons in GeV. For 
P-wave decays the linear momentum dependence is explicitly separated (with p in 
GeV). The amplitudes for a 1.6 GeV glueball should be regarded as our predictions 
for the primitive glueball. The other amplitudes are to be used for calculation of 
the decays of the physical states. The widths (including all partial waves) are listed 
in the final column for a 1.6 GeV glueball, assuming that all resonances are narrow. 
All calculations are for mock meson phase space, except the widths in brackets in 
the last column, which are for relativistic phase space with 70 chosen to agree with 
lattice QCD predictions for glueball decay to tttt, KK and 7777. When both S- and 
D-wave amplitudes are possible the amplitude is the S-wave amplitude. The ratio 
(D-wave amplitude)/^ 2 S-wave amplitude) is -4.2 and -4.1 for decays of a 1.6 GeV 
glueball to pp and lulu, respectively, and -4.0 for decay of a 1.8 GeV glueball to 
K*K*. (irir)s stands for a (hypothetical) narrow {uu + dd) /y/2 resonance /o(600) 
which decays dominantly to titi ||. It may be related to the low mass tail of the 
/o(1370)//o(400 — 1200). Due to the large width and uncertain mass of the physical 
(uu + dd)/ y/2 scalar resonance, the predictions should be viewed as being anywhere 
between zero and the upper limit quoted. 
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